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We compute the time-dependent entanglement entropy of a CFT which starts in 
relatively simple initial states. The initial states are the thermofield double for thermal 
states, dual to eternal black holes, and a particular pure state, dual to a black hole formed 
by gravitational collapse. The entanglement entropy grows linearly in time. This linear 
growth is directly related to the growth of the black hole interior measured along "nice" 
spatial slices. These nice slices probe the spacelike direction in the interior, at a fixed 
special value of the interior time. In the case of a two-dimensional CFT, we match the 
bulk and boundary computations of the entanglement entropy. We briefly discuss the long 
time behavior of various correlators, computed via classical geodesies or surfaces, and point 
out that their exponential decay comes about for similar reasons. We also present the time 
evolution of the wavefunction in the tensor network description. 



1. Introduction 



A black hole interior is, almost by definition, difficult to probe. Perturbatively, an 
outside observer can never receive signals from behind the event horizon. In a unitary 
quantum theory, however, information should not be lost. According to the gauge grav- 
ity duality, all the relevant information is contained in the field theory that lives at the 
boundary. 

In this paper we discuss some observables whose gravity computation involves the 
interior. By 'interior', we mean the region to the future of the event horizon heading 
towards the singularity, not to be confused with the second exterior region of an eternal 
black hole. We consider both the eternal black hole, dual to a thermal state in CFT, and 
a particular black hole microstate that can be formed by gravitational collapse and has 
only one asymptotic region. 

Our main tool is the entanglement entropy, computed holographically as the area of 
an extremal surface in AdS ending on the boundary |1],|2[. In a static situation, extremal 
surfaces do not penetrate the event horizon ||, but if the system is time dependent then 
it is possible to probe the interior. A similar setup has been explored using the Vaidya 
spacetime, which describes a shell of null dust falling into a black hole 0,|],||,|],[7],|| . However 
here we will use only the empty black brane geometry. 

To understand how we introduce time dependence to probe the black hole interior, 
first consider an eternal black hole. It is static under time evolution that runs forward on 
one side of the Penrose diagram, and backwards on the other. However, instead we will 
choose to time-evolve forward in both exterior regions, so the system is time dependent. 
This is a simple model for thermalization in the strongly coupled dual CFT. The two-sided 
setup sounds artificial but in fact it is relevant to a more realistic thermalization process: 
There is a class of black hole microstates that are similar to the eternal black hole outside 
the horizon, but they lack the second asymptotic region. These black holes are dual to 
time-dependent pure states in CFT undergoing thermalization. 

Let us summarize the calculation. The entanglement entropy of a quantum system 
is defined by separating the system into two parts, A and B, on some fixed-time slice. 
The reduced density matrix of subsystem A is computed by tracing out other degrees of 
freedom, pa = Ttb ptotah an d the entanglement entropy is Sa = —TrpA^ogpA- For A 
a spatial region in a two-dimensional CFT, the entanglement entropy in the groundstate 
and a thermal state was computed in ||,|10| using the replica method. 
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Similar techniques were used in [11] to compute the time-dependent entanglement 
entropy of a 2d CFT in a gapped excited state. This is interpreted as modeling the time 
evolution after a global quantum quench, where we prepare the system in the groundstate 
of a gapped Hamiltonian H and suddenly change H — > H, so the system is no longer 
in the groundstate. The system starts with only short-range entanglement. The result of 
TIH is that the entanglement grows linearly, and eventually saturates at the thermal value 



after a time t ~ L/2 where L is the size of region A. That is, the subsystem eventually acts 
thermal due to correlations with the rest of the system. The state that is actually used in 
TUl is given by starting from a boundary state in the CFT and evolving it in Euclidean 



time in order to render it normalizable. We call the resulting state "the B-state" . We will 
point out that the gravity dual of this procedure is very simple. It amounts to cutting the 
usual eternal black hole Penrose diagram in half by adding an end of the world brane in 
the bulk. 

On the gravity side, as mentioned above, the entanglement entropy is related to an 
extremal surface in AdS that ends on the boundary of the region A at the boundary. The 
proposal, made for static spacetimes in |]J and extended to time-dependent backgrounds in 
0, is Sa = Area(7)/(4G/v) where 7 is the extremal surface and Gn is Newton's constant. 
If there are multiple extremal surfaces ending on region A, the rule is to pick the one with 
minimal area. 

We take region A to be half of space, or a finite strip of width L. For the eternal black 
hole, A consists of two identical pieces, one on each side of the black hole. For L ^> (3 
with (3 the inverse Hawking temperature, the minimal area extremal surface at early times 
extends across the black hole from one asymptotic region to the other. As we evolve in time 
it progresses into the interior along special nice slices. After a time t ~ /?, the extremal 
surface gets 'stuck' at a particular time in the interior, but continues to stretch along the 
spacelike direction of the interior. This growth of nice slices along the spacelike direction 
in the interior is directly responsible for linear growth in entanglement entropy, Sa = svt, 
where s is the thermal entropy density and v < 1 measures the speed of entanglement 
growth. 

If A is a finite region of size L, then after t ~ L/2, the minimal area extremal surface 
jumps to a static surface sitting outside the horizon that reproduces the thermal entropy. 
All of these statements also apply to the single-sided black hole microstate. 

In two boundary dimensions, where the CFT calculation can be done explicitly, the 
gravity calculation gives linear growth with v = 1. For the eternal black hole we find an 



2 



exact match to a thermal double CFT. In fact, when the region A is half of space, the 
result is fixed by conformal symmetry. For the single-sided black hole our result agrees 



with the Cardy-Calabrese B- state [11 



The main conclusion is that the growth of the nice slices along the spacelike t-direction 
in the interior is responsible for linear growth in entanglement entropy during thermaliza- 
tion. We also discuss the relevance of this interior region for other observables, including 
correlators of heavy particles, strings, and branes. 

In our setup the entanglement at t = is entirely short rangeB. As the system 
evolves, entanglement spreads out over larger distances. This suggests a picture for the 
wavefunction of the CFT during thermalization, where the ordinary renormalization group 
flow is supplemented by additional structure at scales larger than f3. 

The paper is organized as follows. In section 2, we calculate the entanglement entropy 
from gravity in AdSd+i- In section 3, we specialize to AdSs/CFT 2 , and compute the 
entanglement entropy from both bulk and boundary. In section 4 we consider the region 
in the interior relevant to computing the long time behavior of the expectation values of 
local operators, Wilson loops and higher dimension surface operators. We comment on 
implications for the CFT wavefunction in section 5, and conclude with a discussion of the 
interpretation of the interior in section 6. 



2. The growth of the entanglement entropy 

We first consider black holes that correspond to the holographic duals of thermal field 
theories. We assume that we have one or more infinitely extended spatial directions, and 
consider a metric of the form 



ds' 



-g 2 (p)dt 2 + h 2 (p)dx 2 d _ 1 + dp' 



(2.1) 



We assume there is a horizon at p = and that for small p, 



g(p) = p + o(p 3 ) 



h(p) = C + o(p 2 ) , (C^O) 



(2.2) 



1 This differs slightly from the Vaidya calculations mentioned above, where the initial en- 
tanglement after the quench is spread over different scales. Other types of quenches have also 



been studied holographically in, for example, [12, 13, 14, 1^,1(3, 17, 18, 19]. We can view field theory 
in de-Sitter space also as a particular time dependent state. The growth of the entanglement 



entropy in that case [20] has features similar to those discussed here. 
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Fig. 1: (a) Penrose diagram of the maximally extended black branes we consider. 
We suppress the spatial coordinates along the brane. There are two exterior regions 
Ex and E2, each of which has a boundary. The two copies of the corresponding 
field theories live at these boundaries. There are two interior regions In to the 
future and In' to the past. This can obtained from a Euclidean solution by cutting 
and pasting at a moment of time reflection symmetry, as indicated in (b). This 
procedure also produces the entangled state in the two copies of the field theory. 
It just corresponds to doing Euclidean time evolution over a time j3/2. 



We set the temperature to /9 = 1/T = 2n, unless we explicitly introduce (5. We also assume 
that g, h grow sufficiently fast for large p so that we have a timelike boundary. Simple 
examples of such metrics are planar black branes in AdSd+i where (see appendix A) 

h=- (cosh ^ , g = h tanh ^ . (2.3) 

d V 2 J 2 

The full extended Penrose diagram is shown in fig. 1. (The particular case of AdSs is 
special and will be treated in more detail in the next section). 

The interior region corresponds to p = in and t = tj — itc /2, so that pe l is finite as we 
cross the horizon, and tj is real. The coordinate ti is spacelike in the interior. Note that 
then —iglin) is real and positive in the interior. 

The extended Penrose diagram also contains a second copy of the exterior. The full 



spacetime is dual to the thermofield double pT|]2"2jl . It has two boundaries which correspond 
to two copies of the field theory. These two copies are in an entangled state of the form! 



o J - Jl 



^) = }_^\E n ) 1 \E n ) 2 e-i^ . (2.4) 

n 



The index n is really continuous if we consider field theories in a non-compact space. 
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We will now consider the entanglement entropy of a region "A". Region A consists 
of half of the space of each of the two copies of the thermofield double. We will be 
separating each copy in two halfs at some time tb which is the same on both sides. In 
other words, we take time to run forwards on both copies of the field theory i. Running 
time forwards on both copies we are introducing some time dependence into the problem 
and the entanglement entropy will depend on time. This corresponds to replacing e~^ E " — > 

e -§E n -2iE n t b in 

The conjectured holographic recipe for computing this entropy is the following 
We have to find the area of an extremal codimension two surface that ends at the boundaries 
of the regionB. 

special extremal slice 
extremal surface / 




Fig. 2: (a) Extremal surface that computes the entanglement entropy. For large 
tb, it sits very close to a special critical spacelike surface in the interior. This region 
gives rise to the linear behavior in tb in the entanglement entropy, (b) The same 
for the B-state. Here we can have the surface ending at the end of the world brane. 



In our case, the symmetries of the problem suggest that if region A is the half space 
x\ > 0, then the bulk surface will be extended along X2, ■ • • , %d-i and will sit at X\ = 0. 



3 Of course, running time forwards in one copy and backwards in the other is a symmetry of 
the problem. Our time evolution is not a symmetry of the problem. 

This formula is an unproven conjecture, but it has been shown to be valid in some special 



cases and has passed several checks [23|. 
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Thus finding the surface reduces to finding a function t(p) (or p(t)) which describes how 
the surface moves in the t, p direction. In order to find it, we need to extremize the area 



A = V d . 2 J dt [h(p)] d - 2 y/-g2(p)+fP . 



(2.5) 



We first look for an extremal surface of the the form shown in fig. 2. It has a t j — > —tj 
symmetry under the reflection that exchanges the two sides. So, at tj = we expect that 
p = 0. 

It is easy to find the first integral of the equation of motion of ( |2.5| ). Since ( J2.5| ) is 
time independent, the "energy" is conserved, giving 



2ud-2 



9 2 h 



V~9 2 + P : 



~Woh 



d-2 



(2.6) 



where ho and go are g and h evaluated at the interior point p = ino where p = 0. (Note 
that go is purely imaginary). From here we can find 



tip) 



dp' 



7T 
-I — 

9X ' 1 9 2 h 2 d - 4 



(2.7) 







p complex 






P 


p-0 





Fig. 3: Contour in the p plane that gives t(p). We avoid the pole at p = as 
indicated. 



Note that the integral ( |2.7[ ) has a pole at p = and should be taken along the contour 
indicated in fig. 3. This pole at p = gives the expected behavior t ~ — log p as we cross 

6 



the horizon, implying that pe l remains finite as the trajectory crosses the horizon. The 
value of the area is also simply given as 

/poo h d ~ 2 
dth d ~W-9 2 + P 2 = 2V d - 2 / dp . (2.8) 

Y 1 g 2 h 2d - 4 

The factor of two arises from the symmetry of the configuration, see fig. 2(a). 

In general this integral is finite as long as po = ino is not exactly the point where 
a = (—ig)h d ~ 2 is extremized. Let us denote the point where a is extremized by p m = in m . 
Note that a is zero at the horizon p = and then it starts growing as we go into the 
interior. It typically gets to a maximum value at K m and then decreases again. This can 
be checked explicitly for (|2.3|). 

For finite values of = t(p = oo) we find that Ko < K m . But as becomes large, we 
find Ko K m and that is the reason why tb in ( |2.7p becomes large. For large the minimal 
area extremal surface has two regions. It lies at k ~ n m for a long range of the spacelike 
interior t-direction, and then it crosses the horizon and goes to the boundary. The region 
of the solution where it crosses the horizon and goes near the boundary has a shape that 
depends on the details of the geometry but, for large tb, it gives just a constant to ( |2.8|) . 
On the other hand, the piece lying in the interior gives a large contribution linear in tb to 
the area. This linear contribution is computed by approximating the area by setting p = 
and g — >■ g m , h — > h m in the first equality in (|2T8|). The form of the entanglement entropy 
for large tb is 

S = — — = — — [a m tb + constant] . (2.9) 
Of course A in (|2.8[) is divergent at large p, but the divergence is independent of tb- The 



coefficient a m = \/—g 2 n h d n 2 is the value of a at the extremum, which is a maximum. 



These surfaces were also considered in |M| , here we are expanding on their interpretation. 



It is interesting to express this number in units of the entropy density, s — 
This gives us a measure of the speed of growth of the entropy as 

ds y^glht 2 

— = vV d - 2 s , v = 

9t b h(0) 

For the AdSd+i black brane metrics (|2.3|) we find that 



^ = vV d - 2 s , v = V L/ ^_T • (2-10) 



dK m I d Vd(d-2)* \ 

tan— — = \- — -, v = — t . (2.11) 

2 Vrf-2' ^(rf-l)] 1 -^ 

For d = 2 v = 1, for d = 3 v is slower than the speed of sound and for d > 3 v is larger 
than the speed of sound and always bigger than 1/2. 
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2.1. Other shapes for region A 

Let us now discuss what happens if we take the A region to have different shapes. A 
simple shape is I x R d ~ 2 on each of the two copies of the field theory, where / is an interval 
of length L 3> ft. The minimal area surface consists of two disconnected surfaces, each of 
the form discussed above. Thus we get the same answer, up to a factor of two. This is the 
correct answer for times smaller than L. When the time becomes of order L then there 
is another extremal surface whose area becomes smaller. It consists of two disconnected 
surfaces, each staying outside the horizon, but close to the horizon. For large L their 
contribution goes like the black hole entropy density times the length of the interval. 

More generally, consider a region A given by the union of two identical regions R in 
each of the copies of the field theory. We take R to be very big compared to f3 (and assume 
that its boundary is not varying rapidly on distances smaller than its size). Then we get 
the same result as in ( |2.9| ) with Va-2 — > Bd-2 where Bd-2 is the area of the boundary of 
region R. 

Notice that the thermofield double ( |2.4| ) (at t = 0) is behaving as a state with a mass 
gap from the point of view of the two copies of the field theory. Here what we mean is 
that all equal time correlation functions decay exponentially as we separate them in space. 
Also the entanglement entropy for the regions in the above paragraph goes like the area 
of their boundaries. Of course the system is not gapped as a lorentzian solution, where we 
can have very low energy excitations by considering particles very near the horizon. 

When tb becomes comparable to the size of the region R, we can have other surfaces 
which can give a smaller area. These are surfaces that stay on each of the exterior regions. 
The surfaces are disconnected, they live purely at tb and do not cross the horizon. However, 
they asymptote to the horizon and, for large regions, they give an entropy of the form 
S = 2Vrs, where s is the entropy density of the black brane and Vr is the volume of the 
region R. 

2.2. Black holes formed from pure states 
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Fig. 4: (a) Penrose diagram for the B-state. We have an end of the world brane or 
orbifold plane along the dotted line. This breaks the time translation symmetry. 
The corresponding state can be produced by cutting a Euclidean solution at a 
moment of time reflection symmetry, as indicated in (b). Note that in the Euclidean 
solution the end of the world brane reaches the boundary. Thus the Euclidean field 
theory has a boundary. The state at t = corresponds to performing Euclidean 
time evolution on the conformally invariant boundary state by an amount /3/4 of 
Euclidean time. 



The behavior of the entanglement entropy under black hole formation was studied in 
numerous publications which found features which are the same as what we discuss here 
[|^,^|,|^,|6|,[T|,|^, p73] , pr4| , p7^ , p~6] , p7 : 7l , p~8[1 . Several of these publications focused on the case of Vaidya 
metrics. Here we would like to discuss a very simple black hole formation situation which 
is related to the above discussion. In addition, it is the precise holographic dual of the 



idealized quench configurations discussed in ill 



The solutions are simply the eternal black hole solutions with an end of the world 
brane that cuts them in half. This end of the world brane can have different physical 
origins. It could be a true end of the world brane in an M-theory solution, or it could arise 
from various types of orbifolds. Some examples are given in ||22||i. 



5 An example of an end of world brane is the end of the world brane of M-theory. We can 
get others from type F theory or various orbifolds. We do not want to consider here the "geon" 



configuration discussed in [25], because those solutions involve a boundary condition which is 
not local in the spatial coordinates of the CFT (even though they are perfectly well defined and 
reasonable objects to study). 
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This is a state whose description is almost as simple as the one above. Namely, we 
consider a field theory with a boundary condition. We then evolve in euclidean time for a 
time /3/4, 

\y} = e-* H \B}, (2.12) 

where \B) is a boundary state in the conformal field theory. We can start with this state 
at t = and evolve it in Lorentzian signature. 

The state at t = has a mass gap of the order of the temperaturj. This is a common 
feature of (stable) AdS black holes. 

The holographic computation of the entanglement entropy is essentially identical to 
the one we had above. The extremal surface used to compute the entropy can end on the 
end of the world brane. This end of the world brane sits at tj = 0. So the solutions are 
the same as the ones we considered above. And the entanglement entropy is essentially 
the same, up to a factor of two, since we now consider only one of the sides. 

We see that as time evolves the entanglement entropy grows. This growth of the 
entanglement entropy can be understood from two points of view. From the field theory 
point of view, the local entanglement that existed at t = is spread out by the time evo- 
lution. Thus, for large times the entanglement is highly non-local along the non-compact 
directions of the field theory and it is spread over a region of size t}>- 

In the bulk solution the entanglement is given by a surface that starts from the 
boundary and ends on the end of the world brane. See fig. 2(b). The linear growth in % 
is due to the fact that the extremal surface is growing along the ti direction inside the 
horizon. The position along the timelike direction in the interior is barely changing. 

The shape of the extremal surfaces is essentially identical to the shape of the so called 
"nice slices" that are usually introduced in discussions of black hole information loss [2t| . 



The growth of the entanglement entropy is due to the growth of these nice slices. These 
slices grow because we keep adding space in the interior, making the range of tj bigger. 



6 Again, this is a mass gap in the sense that the equal time correlation functions decay 
exponentially in space. Also, the entanglement entropy of a region goes like the area of the 
region. Of course, the lorentzian spacetime can support very low energy excitations which are 
obtained by placing particles very near the horizon. 
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Fig. 5: Euclidean cylinder with twist operator insertions used to compute the 
entanglement entropy Sa in a finite-temperature CFT, where region A is (a) the 
half-line, and (b) a finite interval. 



3. BTZ black strings and two dimensional CFT 

Here we consider a two dimensional CFT on an infinite spatial line. In this case the en- 
tanglement entropy can be computed in CFT using the replica technique ]^|| 0^ , EM|23 ] 
and can be matched quantitatively to the gravity formula. We will also describe the gravity 
calculation in the previous section from an alternative viewpoint, using the fact that the 
'planar' black hole in AdSs is just ordinary Ad S3 written in different coordinates. 

We consider four scenarios, both in gravity and in CFT: the subregion is either a half- 
line or a finite interval, and the CFT is either in a thermal state or in the pure B-state. 

3.1. CFT at finite temperature: half-line 

First we consider a finite-temperature CFT and choose region A to consist of the half- 
line x > 0, both in the CFT and in its thermal double. In Euclidean signature, insertions 
in the CFT are at Im z = and insertions in the thermal double are placed at Im z = i/3/2. 
At time t = 0, the n th power of the reduced density matrix, p 1 ^, is a Euclidean path integral 
on an n-sheeted cylinder with z ~ z + i/3, shown in fig. 5(a). This is given by a two-point 
function of twist fields inserted at the branch points. The answer at time t is computed 
by analytically continuing the positions of the twist fields to Lorentzian time, 

I n = trp"X = (^ + (z 1 ,z 1 )^~(z 2 ,z 2 )} 

(3-1) 

z x = zi = , z 2 = -z 2 =2t + i(3/2 . 

Note z\ 7^ z 2 because of the continuation to Lorentzian t. We have chosen time to run 
forward in both copies of the CFT (giving the factor of 2), so it is not a symmetry of the 
thermal state. The twist fields have dimension 

*. = £(„-!), (3.2) 
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where c is the central charge. 

Mapping to the plane w = e 27rz ^ gives the two-point function 

(3 , {2nty~ 2A 
-cosh — — 



27T6 U i 



(3.3) 



where the short-distance cutoff e comes from regulating the twist fields. Therefore the 
entanglement entropy Sa = —d n I n \ n =i is 

S A = ^ log (^2 cosh *j\ + 2S div , S div = I log . (3.4) 

We see that for t > j3, the entanglement entropy grows linearly in time. 

3.2. CFT at finite temperature: finite interval 

Now take region A to consist of two identical intervals of length L, one in the CFT and 



one in the thermal double. The case t — was discussed in |30| and is drawn in fig. 5(b). 
At finite time, the reduced density matrix is given by a 4-point function of twist fields, 

I n =($ + (z 1 ,Z 1 )$-(z2,Z2)$ + (z3,Z 3 )$~(z 4 ,Z 4 )} , 

zi = zi = , z 2 = z 2 = L, (3.5) 
z 3 = 2t + L+Z/3/2, z 3 = -2t + L - 0/2, z 4 = -z 4 = 2t + 0/2 . 

Mapping to the plane, this becomes 

o \ -4A„ / 7rL\ _4A " a 

2 sinh — — — ) (^cross^-cross) G n (x cross , x cross ) (3-6) 



2nJ \ (3 

where x croS s is the conformal cross-ratio 

x _ (wi - w 2 )(w 3 - w 4 ) 2sinh 2 (7rL//?) 

{w\ — w 3 ){w2 — w 4 ) cosh(27rL//3) + cosh(47rt//?) 

and 

G n {w,w) = ($+(0)$"(«;,«;)$+(l)$-(oo)) p;ane . (3.8) 

We also see that a; cr0 ss = Across- At this point we have mapped the problem to the 4-point 
function of twist correlators on the plane, which is the groundstate entanglement entropy 
of two disjoint segments in a single CFT. The extra factors in ( |3.6| ) come from the fact that 
the two intervals are boosted. The problem on the plane in theories with a gravity dual 
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was studied in detail by Headrick [31], building on previous work on multiple intervals 
in |3^,[J3]]. Therefore our calculations are just a small variation on |3l[ with a different 
interpretation because now the cross-ratio depends on time. 

At early times /5 C f < 1/2, we see that we are doing the OPE around x cross ~ 1 
and the correlator factorizes, 

I n « ($ + (z 1 ,Z 1 )$-(z4,Z4)}($ + (z2,Z 2 )$~(z 3 ,Z 3 )} 

a \"4A„ (3.9) 
' \ e -8irA n t/p _ v ' 



4ive 



Because of the exponential dependence on t,L, this formula is actually valid until t is 
within {3 of L/2. We see that we get twice the answer for the half line in (3TT), for any 



CFT. At late times t ^> L/2, the correlator factorizes in the other channel x cross ~ 0, 

In^[£^j ^" e" 47rA " L ^. (3.10) 
From these expressions we find the entanglement entropy 

[ 4S div + ^ t < L/2 

Sa= { 2ncL ' (3 - n) 

[4, div + ^ t>L/2 

with ,5 div in 

That is, the entanglement entropy grows linearly until it saturates at the thermal 
value. Generically, the transition between the two regimes is smoothed out over a time of 
order (3. In theories with gravity duals we get a sharp transition, discussed below. 

3.3. CFT in the pure B-state 

Now consider a CFT in a pure state 

|^ = e -Pm\B) , (3.12) 

where \B) is a boundary state. Correlation functions in the state |^) are computed by a 
Euclidean path integral on a strip, 

< Im z < - , (3.13) 
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with the boundary condition B imposed at the top and bottom of the strip. Therefore the 
entanglement entropy of a half-line at time t is related to the one-point function on the 
strip, 

I n = ($+(*!, z x )) , z 1 = -z 1 =t + i/3/4 . (3.14) 

By mapping to the upper-half plane w = e 27rz ^ where the conformally invariant one-point 
function is (w — w)~ An , we find 

Therefore the entanglement entropy of the half line, Sa = —9 n I n \ n -i, is 

Sa = % log ( 2 cosh — — ) + S'div • (3.16) 
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This is half of the answer obtained in ( |3.4|) for the half-line in a thermal state. The reason 
for this is clear: in the pure B-state, when we compute the correlator on the upper-half 
plane, there are image points in the lower-half plane. The image points play the same role 
that insertions in the thermal double played in the previous calculation. The same is true 



for the finite interval (see [|lT| for details of the computation) in the two OPE limits that 



we considered. In summary, we have in both cases 

gpure B— state £ ^thermal double yj^ 

for the linear rise and, for the interval, the late-time constant behavior. 
3. 4- BTZ black string 



The BTZ black hole is a quotient of AdSs [p4]| . The results of section || apply to 



black holes with an infinite spatial direction, so let us unwrap the BTZ angle, giving back 
ordinary AdS^. Obviously now there is only one boundary. When we unwrap the quotient, 
the two boundaries of the BTZ Penrose diagram become the two wedges of Rindler space in 
the CFT |35|] (see |3(J for a more recent discussion) . This maps the problem of computing 



extremal surfaces in section |2] to a very simple one, so we can find the extremal surfaces 
explicitly and see exactly how they pass through the bulk. 

First we will describe the various regions of the black string. It is useful to relate each 
coordinate patch to Poincare coordinates, 

1 

7 2 



ds 2 = -z{-dxk + dxj + dz 2 ) , (3.18) 
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Fig. 6: Left: Bulk and boundary regions of the BTZ black string. The boundary 
regions I and III are the two Rindler wedges, and region II is the Milne patch. The 
exterior of the black string ends on region I and extends perpendicularly into the 
bulk. The shaded region in the bulk is the interior, which reaches the boundary 
only along the light cone x\ — x\ = 0. The future bulk region lies to the future of 
the shaded region and its boundary is the Milne patch. Right: Endpoints of the 
entanglement region. 



which ends on the Minkowski diamond of the CFT, pictured in fig. 6. The exterior metric 
of the black string is 

Exterior : ds 2 = — sinh 2 pdt 2 + cosh 2 pdx 2 + dp 2 . (3.19) 

These coordinates cover a portion of the Poincare patch. Near the boundary, they are 
related by 

x x ± x ~ e~ x±t , - « -e p ~ x , (3.20) 

% 2 

so we see that t, x cover one Rindler wedge of the Minkowski diamond. This is region I in 
fig. 6. We can reach the other Rindler wedge (region III) by continuing t — > t + in, which is 
the same continuation that would take us from one side of the black brane Penrose diagram 
to the other. 

The black string metric (|3.19| ) has a horizon at p = 0. To go behind the horizon we 
continue p = ia, t = t — in/2, and find the interior metric 

Interior : ds 2 = sin 2 adt 2 + cos 2 adx 2 — da 2 . (3-21) 
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These coordinates cover the shaded interior region in fig. 6, which meets the boundary 
only along the light cone x\ — x\ = 0. Note that the interior of the black string does not 
correspond to the future diamond of the Milne parameterization of flat space, region II. 
That bulk region, to the future of the shaded portion in fig. 6, can be reached by setting 
a = it/2 — ip,x = x — z7r/2, so 

Future : ds 2 = cosh 2 pdt 2 — sinh 2 pdx 2 + dp 2 . (3.22) 

The three bulk regions that we have just described — exterior, interior, and future — are 
related to Poincare coordinates in appendix B. The boundaries of the different regions are 

Exterior : X 2 , — x\ < 

Interior : < x\ — x\ < z 2 (3.23) 
Future : z 2 <x\ — x\ . 

There are actually two interior patches, inside the past and future horizons; the coordinates 
( |3.21[) cover the patch with xq > 0. 

3.5. Holographic entanglement: half-line 

Now, consider a subregion A that consists of the half-space x > in each of the 
two Rindler regions, at time tb- We will apply the holographic entanglement formula to 
compute Sa- In Poincare coordinates (a?o,xi), the endpoints of region A in the left and 
right Rindler wedges are 

Pi = (sinh tb, — cosh tb) , Pi = (sinh tb, cosh tb) ■ (3.24) 

The two points are at equal times and separated by Axi = 2 cosh tb, see fig. 6. Note that 
the computation of the entanglement in a CFT and its thermal double has been mapped 
to the problem of a single interval in a single CFT. The geodesic that joins Pi and P% is 
just a semicircle, given by the equations 

z 2 + x\ = cosh 2 tb, xq = sinht;, . (3.25) 

We will follow this semicircle through the various regions of the black string. In the 
coordinates (|3.19| ), these equations imply x = and 



sinhttanhp = sinhtb . (3.26) 



16 



We see that as p decreases, t increases and goes to infinity at the horizon. In the interior 
( p. 2 1|) , the semicircle obeys 

cosh t tan a = sinh tb . (3.27) 

As a increases from zero, t decreases until t = where it meets the curve coming from 
the other Rindler wedge. If tb 3> 1, then a will be very close to a m = 7r/2 for most of the 
range t < t b . If we had performed the BTZ quotient, this special slice at a m would be the 
singularity of the BTZ black hole. 

We compute the length of the semicircle in Poincare coordinates, 

^coshtf, t 

A 12 = 2cosht b ^^=^^= . (3.28) 



o 



z 



a/ cosh 2 tb — z 2 



The divergence near the boundary is regulated by placing a cutoff at p = p m ax- Defining 
e = e~ Pmax and using ( |3.20| ), the integral is cut off at z = 2e. The final answer for the 
entanglement entropy is then 

S, = ilog(^). (3.29) 
Reintroducing the temperature f3 = 2n and Newton's constant Gn — this becomes 



Sa = ~ log (^2 cosh y 1 + 2S div , (3.30) 



in agreement with the CFT calculation ( |3.4| ) 



3.6. Holographic entanglement: finite interval 

Now let us choose region A to be an interval of length L on each side of the Rindler 
diagram. The endpoints are shown in fig. 6. In Minkowski coordinates, we see that this is 
just the calculation of the entanglement entropy for two boosted intervals. 

There are two ways to draw the bulk region ending on region A, as shown in fig. 7. 
(The figure is projected onto the Xi,z plane; since the intervals are boosted, it is not a 
single time-slice.) Assuming tb f3, the minimal length surface at early times consists of 
two geodesies that cross from one Rindler wedge to the other: one from Pi to P 2 that was 
described above, and one joining the points 

P 3 =e" L Pi, P 4 = e~ L P 2 . (3.31) 
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Fig. 7: Geodesies that connect the endpoints of the entanglement region, pro- 
jected onto the 2:1,2: plane. There are two choices, shown as dashed or dotted 
lines. Dashed geodesies go through the interior and dominate at early times. Dot- 
ted geodesies stay within each exterior region and dominate at late times. 



These are the dashed lines in fig. 7. The two geodesies have the same regulated length, 
because the L-dependence disappears when we use the appropriate cutoff p = Pmax, z = 
2ee L . Therefore we just multiply ( |3.2y| ) by two, and find 

S { a } = ^h + 4S div (3.32) 



in agreement with the CFT result ( |3.11[ ) at early times. 

The other configuration of geodesies ending on region A (drawn as dotted lines in 
fig. 7) is to connect Pi to P3 and P2 to P4. These geodesies do not pass through the 
interior. Instead they skim the outside of the horizon, so we expect an answer related to 
the horizon area, i.e., the thermal entropy. It is easiest to work in exterior coordinates 
( |3.19|) . The symmetries of this configuration allow us to shift x — » x — L/2 so that the 
geodesies lie in equal-t slices, then set t = because the metric is static. Now we compute 
the distance between 

^ = (e- L / 2 ,0), Pl = (e L / 2 ,0). (3.33) 
The geodesic is again a semicircle in Poincare coordinates and its length is 

A 24 = cosh (%) [ - . 1 , (3.34) 

W J z ^cosh 2 (L/2)-z 2 

where the integral is over two segments: one from z m i n = 2ee L / 2 to z max = cosh(L/2), and 

one from z m i n = 2ee _L / 2 to the same z max . Reintroducing (3,G^ as above and assuming 

(2) 

L 3> {3, the entanglement entropy S A = (A13 + A2a)/A is 

Sf = + 4S div . (3.35) 
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Comparing to the other configuration ( |3.32|) , we see that indeed the minimal area surface is 
the connected diagram at early times and the disconnected diagram at late times. Putting 
everything together, we have linear growth in the entanglement entropy given by S A ^ until 

(2) 

tb = L/2, followed by a sharp transition to the thermal value S A . This agrees with the 
CFT result (|3.11| ), but now the transition is sharp. 

The fact that there is a sharp transition rather than a smooth cross-over in CFTs with 
a gravity dual was understood in |3lj]. The configuration in [31] is two disjoint intervals 



[0, x cross ] U [1, oo] at fixed time, and the sharp transition appears at x cross = \. It is argued 
that this transition is related to the Hawking-Page phase transition between thermal AdS 
and the BTZ black hole (and higher-genus analogues of this transition). To see this, note 
that in the case n = 2, the replica CFT lives on a torus, so sharp transitions in the torus 
partition function can be expected to lead to sharp transitions in the entanglement entropy. 
(See |JT| for details.) In our case the intervals are boosted, and the cross-ratio in ( |3.7| ) is a 
function of time. The transition point x cro ss = \ corresponds to t = L/2, so this becomes 
the transition between Sa ~ t and Sa ~ constant. 

3. 7. Holographic entropy of the pure B-state 

In the case of the single sided black string described in section |2.2| , the answer for the 
entanglement entropy ( p.30|) , ( |3.32| ), or ( |3.35| ) is simply divided by two. The holographic 



dual of the single sided black string is time evolution in a pure state [*) = e- pH / A \B), 
which is precisely the Cardy-Calabrese calculation described in section There we 

also found the entanglement entropy was related to the thermal calculation by a factor 
of two. Thus the single-sided black string also agrees with CFT, for both the half-line 
and the interval, and provides a precise holographic dual of the Cardy-Calabrese B-state 
computations. 



4. A comment on the long time expansion of correlators 

In the previous sections we pointed out that the entanglement entropy had a piece 
linear in time which was produced by a surface living in the interior at some particular 
interior-time p = iK m . Its extension in the spacelike tj-direction in the interior gives rise 
to the piece linear in time. In this section we ask whether the same behavior is present 
for other observables. Indeed, we will find a similar behavior for correlators of operators 
whose duals are heavy particles, strings or branes. 
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First we note that in the case of the BTZ black hole in AdSs, the minimal length 
computation we did for the entanglement entropy is formally identical to the one we would 
do for computing the correlator of two operators, one on each boundary, in the regime 
where we use the classical geodesic approximation in the bulk (A ^> 1). The fact that 
the action is linear in time leads to the exponential in time behavior of the correlator 
(Oi(tb)02(tb)) ~ e ~ 2tbmRAdS 3 . This decay rate is simply the lowest quasinormal mode 
associated to the corresponding bulk field. 

In fact, in any dimension, in the large mass regime, the lowest quasinormal mode is 
related to a geodesic that sits at a particular (possibly complex) value of p in ( |3.19p , and is 
extended along the time direction. The value of p is determined by extremizing the action 
for a geodesic 

iS = - [ dt^-g 2 + p 2 (4.1) 



after setting p = 0. Thus we extremize y—g 2 - In general the extrema lie at complex 
values of p. For example for the black brane metric ( |2.3|) we find that 



cosh(dp c ) = -(d - 1) , y^ = ±e ±i ?c d (4.2) 
where ca is a positive function of d @. These values agree with the quasinormal freque 



lcncics 



discussed in great detail by Liu and Festuccia in |p7j . In this case we get a complex value 
of p c so we cannot say that the geodesic is in the interior. What we wanted to highlight 
here is the simple fact that the geodesic is extended along the t-direction. The saddle 
points in ( ^1.2|) are complex because the function \J —g % does not have a maximum in the 
interior of the black hole, it grows from the horizon to the singularity. 

If one considers particles with imaginary momentum, p = ikm, with real fc, then for 
< k < 1 there are quasinormal modes which correspond to particles sitting in the interior 
at real values of the interior time coordinate k (p = in). As k — > 1 these geodesies move 



closer to the horizon ||37|| . 



Note that the computation of quasinormal modes can be done by solving the wave 
equation completely outside the horizon. Nevertheless the actual answer is well approxi- 



mated by geodesies that lie inside. This was explored in detail in [38,39]. In particular the 



7 Cd = 2 (hp) d Y d^2' which differs by a simple factor from eqn (5.3) in [37], due to a 
different choice of the temperature. In evaluating v/ — g 2 we need to take a dth root, which we 



did following the prescription in [37[. 
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solution of the problem outside dictates which geodesies should be included. On the other 
hand, in the entanglement entropy computation that we did above, we included only the 
real geodesic but no other possible complex minimal surfaces. Presumably a more detailed 
analysis of possible geodesies would show that the real ones are the dominant ones. This 
is a gap in our reasoning that we leave to the future. We have also ignored the possibility 
that the surfaces could hit the singularity. 

We now consider expectation values of Wilson loops and higher dimension surface 
defect operators. More precisely, we imagine correlators of straight surface operators, 
extended along the spatial directions, and localized in the time direction. We consider two 
such operators, one in each of the two copies of the thermofield double at equal times tb- 
Then to extract the large tj, behavior, we need to consider a surface extended along the 
spatial directions and along the (spatial) ^-direction of the interior. To find the interior 
position in the timelike direction we need to minimize 

zS = -Ty^{pj[h(p)} p (4.3) 

where p is the dimensionality of the defect operator, p = 1 for Wilson loops, p = for 
ordinary local operators, etc. We see that the computation of the entanglement entropy is 
similar to a computation of a defect operator correlator with p = d — 2. Minimizing ( |4.3| ) 

we get 

cos dn c = 1 — - — - — r . (4-4) 

The solutions are real as long as p > (for the equal sign we have surfaces that approach 
the singularity). Note that for d = 4 and p = 1 we get a surface that is very close to the 
singularity. 

The expectation value of the defect operator then goes like 

(W) ~ e - 2tb( -- lS ^ (4.5) 

where S c is the value of ( |4.3| ) at (|4.4| ). A similar expression gives the decay of the expec- 
tation value of a defect operator for a black hole that forms from collapse, except that 
2tb — >• tf,- 
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5. Tensor networks and time evolution 



In this section we suggest a tensor network representation for the wavefunction of the 
thermofield double or the pure B-state which is inspired by the above computations. A 
tensor network is a representation of the wavefunction for a discrete quantum system that 
captures the relevant subspace in the Hilbert space where the wavefunction in question 



lives. This has been discussed both for massive theories and conformal field theories [50]. 
It incorporates efficiently the information about the renormalization group. 



Fig. 8: Tensor network representation of the wavefuction of a gapped two di- 
mensional system. The vertical open lines represent each of the spins of a one 
dimensional chain. The intersection points with the horizontal line represent ten- 
sors. The horizontal line gives the pattern of index contractions of these tensors. 
For each value of the spin (the index of the vertical line), the tensor has two indices, 
the two horizontal lines coming out of the vertex. These are contracted with the 
neighboring tensors. 



The simplest case is a spin system in one space dimension in a state with a mass gap. 
The tensor network representation of the wavefunction (called a matrix product state in 
this case |[41|| ) has the form 

*{8 U ---,s L )=Tr[T Sl T aa ---T aL ] , (5.1) 

where T s . = (T s .)f is a D x D matrix and Si = ±1 are the values of the spin at each 
site. See figure fig. 8. This can be viewed as a particular ansatz for the wavefunction that 
parametrizes the subspace of the Hilbert space where we expect the groundstate wavefunc- 
tion to live. It is a space of dimension 2LD 2 which is smaller than the 2 L -dimensional total 
Hilbert space of the system. We need log-D comparable to the entanglement entropy of 
half the system in order to get a reasonable answer. This allows for a drastic reduction in 
the size of the relevant Hilbert space and has proved very useful for the efficient numerical 



computation of the groundstate wavefunction in Id gapped systems [|2| (see []4j| for a 
recent review). 
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Fig. 9: Pattern of index contractions for the tensor network representation for 
the ground state of a massive theory in two spatial dimensions. 




Fig. 10: Tensor network representation for a conformal theory in the IR. At 
each intersection point we have a tensor with five indices and the pattern of index 
contractions that is indicated by the diagram. 



One can write similar ansatze for gapped systems in d spatial dimensions. In that 
case the tensors T will have 2 d indices and their contraction gives a structure similar to a 



d dimensional lattice, see fig. 9, [fH 



For systems that are conformal in the IR there is also a convenient representation 



which is called the MERA network |45[ . It is a particular tensor network with contractions 
as indicated in fig. 10. For more details see | 4"T]|j . This is a representation of the wavefunction 
which incorporates the renormalization group. More precisely, it can realize explicitly a 
discrete subgroup of the renormalization group. As pointed out in [46], we can think of 
these tensors as the wavefunctions of regions with size of order the AdS radiusi. These 



This represents locality well for distances large compared to the AdS radius, but not for 
smaller distances. After all, such small-scale locality is not valid for an arbitrary CFT. 
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tensor networks nicely capture some features of the entanglement entropy, as discussed 
in detail in |46| . |47| , ^8| ,49|. For the Id chain (fig. 8), the entanglement entropy of half the 
system is bounded above by log D, where D is the dimension of the link that must be cut 
to separate the system into two parts. In a more elaborate tensor network, if we must cut 
£ links to separate subsystem A from the rest of the tensor network, then 



S A <£\ogD 



(5.2) 



The MERA network is a discretized version of hyperbolic space, and we can separate 
system A by cutting the links on a minimal-size surface extending into the tensor network. 
Therefore, if the tensors at each site are sufficiently random for the bound ( |5.2|) to be 
saturated, then this resembles the Ryu-Takayanagi formula for holographic entanglement 



entropy [[46 . 




XX XX XX XX XX XX XX XX XX XX XX XX XX 

Fig. 11: Pattern of index contractions for the tensor network representation for 
the thermofield double of a scale invariant theory. The scale invariance is broken 
by the temperature. At scales of order the temperature we have a pattern similar 
to the one we had for the massive theory in fig. 8. 



The thermofield double or the S-state are gapped states with a scale invariant UV 
region. Thus the tensor network representation of those states are expected to be as in 
fig. 11. We essentially have the MERA ansatz up to the scale of the temperature and then 
we have a gapped state. This is the tensor network representation of the state at t = 0. We 
can now ask, what is the representation of the state at larger times t. Our entanglement 
computation leads to a natural proposal. We need to use a tensor network that grows as 
t grows. We simply add more layers in the middle region. We expect the number of these 
layers to be of order t/(3. The resulting tensor network in this middle region looks similar 
to the one we would obtain for a gapped system in d spatial dimensions. The structure of 
the network is similar to the structure of the nice slices. The gapped region is the interior 
of the horizon. The fact that it is gapped is associated to the fact that all correlators decay 
exponentially in time. This representation of the wavefunction of such systems does not 
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Fig. 12: Pattern of index contractions for the tensor network representation 
for the thermofield double of a scale invariant theory after time evolution. As we 
evolve we add of order t//3 intermediate layers. This picture is suggested by the 
geometry of the interior and general properties of time evolution. The red dotted 
line cuts the system into halves, and the entanglement entropy is bounded by the 
number of black lines that it crosses. 

appear to have been used in the tensor network literature, but the black hole suggests that 
it should be useful. 

At some level the representation in terms of tensors as in fig. 12 is obvious, since 
we can view each additional row that we add, as we time evolve, as the action of the 
Hamiltonian on the block spin variables at scale (3. The quick thermalization is implying 
that the resulting state is similar to that of a gapped system, at least for simple local 
operators (the full theory is still unitary). 

We know that for long times we cannot have a perfect decay for all correlators, since 
that would imply information loss. 

Notice that this tensor network idea is not good enough, so far, to resolve distances 
smaller than the AdS radius. The size of the time-like p-direction in the interior is of the 
order of the AdS size and it is not visible in this ansatz for the wavefunction. On the 
other hand, the spacelike t-direction in the interior has a very large size, and is visible. 
Unfortunately since we cannot resolve the timelike p-direction, we cannot say whether we 
are talking about the interior or the region near but outside the horizon. 
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6. Discussion 



Note the peculiar nature of the entanglement at t = 0. At t = the entanglement 
is relatively short range in the non-compact spatial directions. Degrees of freedom are 
entangled within a distance of order (3. This is true both for the two sided and the one 
sided configurations. For this reason the entanglement entropy obeys an area law, similar 
to the one observed in theories with a mass gap. In fact, it is well known that thermal field 
theories behave as systems with a mass gap, when we consider equal time observables. 



Lorentzian copy 1 
Euclidean (3/2 { 
Lorentzian copy 2 

t 













\ / i\ 
/\ 1 \. 
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Fig. 13: Quasiparticle picture of growing entanglement in the thermal double 
CFT. Particles connected by a dashed line in the Euclidean region are entangled. 
One particle of the entangled pair is in each copy. At time = these correlations 
are localized around x ~ /3/2. As time passes the particles move in the two copies 
of the CFT. Particles that contribute to the entanglement are those where one 
member of the pair is in region A and the other is outside. For example the pair 
with the dashed blue line does not contributed because both particles are in region 
A. 



A 

Lorentzian 

Euclidean |3/4 { 

Boundary condition B 

Fig. 14: Quasiparticle picture of growing entanglement in the pure B-state. 
Euclidean time evolution creates entangled particles separated by x ~ /3/2, heading 
in opposite directions. As time evolves one member of the pair can leave the region 
A. This will increase the entanglement entropy. Here the pair with the dotted blue 
line does not contribute to the entanglement, but the rest do contribute. 
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As time evolves this short range entanglement spreads out along the spatial directions. 
This is what gives rise to the linear growth of the entanglement entropy in the boundary 



theory. In fact, we have the picture in fig. 13 and fig. 14, discussed in ||11|| . This picture 
is rather obvious from the boundary theory point of view. It is interesting that in the 
gravity description the same linear growth comes from the linear growth of the spacelike 
t-direction in the interior region of the black brane geometry. This linear growth in time 
and the fact that it comes from the interior region was pointed out in previous work 
@,|,|J|,^,|3],mi|,|6],0,|l|]. We hope to have clarified and isolated the origin of this 
linear behavior as arising from a very special spatial slice in the interior. It is the slice that 
maximizes the area for a codimension two surface. It is also noteworthy that fig. 13 and 
fig. 12 look similar. In fact if we look at fig. 13 and we view the lines as representing index 
contractions and the intersections as representing tensors, then it looks rather similar to 
fig. 12. 

In the boundary theory, when we disturb the fluid with the insertion of a local operator, 
we create a local disturbance that is then carried to longer and longer distances. The 
above results for entanglement entropy suggest that when the excitations have reached a 
size bigger than of order /3, then we should think of them as living in the interior. As 
time progresses these excitations live further and further away in the interior, as measured 
along a spacelike "nice" slice. Since its information is spread out over a long distance it 
seems reasonable that it would be hard to extract via simple local operators. Of course, 
this is not explaining why this excitation cannot get out of the black hole. To understand 
bulk causality, we need to understand bulk locality at distances less than the AdS radius. 

Throughout this paper we have assumed that we have a non-compact direction. If we 
have the QFT on a space of size L, then we can follow the growth of the entanglement 
entropy in this fashion up to a time of order L. On the other hand, the interior seems 
to continue making sense for times which are much longer. It is probable that the "long 



string" phenomenon []50| is relevant for enabling the spread of the entanglement. In other 
words, it is likely that the interior is associated to the spread of entanglement, but that such 
a spread is occurring among the matrix degrees of freedom rather than along the physical 
spatial dimensions. An extreme example could be the black hole in the gravity dual of the 
DO brane quantum mechanics (or any other gravity dual of a quantum mechanics theory). 
In this case, there are no extra spatial directions and the holographic direction, as well as 
the interior of the black hole, has to come solely from the structure of the matrix degrees 
of freedom. 



27 



Recently it was suggested that black holes which are maximally entangled with radia- 



such that its microstates, given by the CFT on one side, are maximally entangled with a 
second system, the second copy of the CFT. The two CFT's do not talk to each other and 
we can view the second CFT as living in the memory of the quantum computer that has 
processed the outgoing radiationi. In this particular case, we get that the horizon is indeed 
smooth. Of course, this state is very special. One can wonder what happens if we applied 
a generic unitary transformation to the Hilbert space of the second side. Could that create 
a firewall? One unitary transformation that we can apply is simply time evolution. We see 
that time evolution does not tend to destroy the horizon. On the contrary, it is making 
the interior region even bigger, in the sense described above. Namely, time evolution gives 
rise to the spacelike t-direction in the interior region, labeled In in figure fig. 4. However, 
this does not address the general case. 
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Appendix A. Change of coordinates for the black brane 

The usual coordinates for a black brane for AdSd+i have the form 



tion do not have smooth horizons ]51| . The eternal black hole is an example of a black hole 



ds 2 = -(l-z d )dt 2 + 



dz 2 



+ dx 2 



(A.l) 



z 



Here the temperature is not l/(27r), but we will fix that later. We now write 



dp 



dz 



tanh — 




d 



— z 



(A.2) 



z 



Then h and g are as given in ( |2.3| ) , where we rescaled t and x to introduce the factor of | 
to fix the temperature to l/(2ir). 



It has been claimed that it is impossible to construct such a computer [52]. 
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Appendix B. Coordinates of the BTZ string 

In this appendix we relate the various coordinate patches of the BTZ string to Poincare 
coordinates via the embedding space 

-Yl x - y 2 + if + Yi = -i . (b.i) 

The Poincare coordinates ( |3.18| ) are 

Y_ 1 = ±[l + (z*-xt + x 2 1 )l Y 2 = ±[l-(z*-xt + xl)}, y = -, yi = -.(B.2) 

zz Zz z z 

The exterior coordinates ( |3.19| ) are 

Y_i = cosh x cosh p, Y2 = — sinh a; coshp, Yq = sinh £ sinh p, Y\ = cosh £ sinh p , 

(B.3) 

which covers — 1q 2 + ^i 2 — 0? ~^-\ + ^2 2 — — 1- The interior coordinates ( 3.21 ) are 



Y"_i = cosh x cos a, Y2 = — sinh x cos a, Y"o = cosh £ sin a, Y\ = sinh £ sin a . (B.4) 

The interior covers -Y 2 + Y? < 0, -Y\ + Yg < 0, or -1 < -Y 2 + Y? < 0. Finally the 
future region ( |3.22| ) is related to the embedding space by 



Y"_i = sinh x sinh p, Y2 = — cosh x sinh p, Yq = cosh £ cosh p, Y\ = sinh £ cosh p . 

(B.5) 

This region covers -Y 2 + Y 2 < -1, —Y\ + Yg > 0. 
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